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Abstract 

On the basis of a new ab initio, all-electron response scheme, formulated within time- 
dependent density-functional theory, we solve the puzzle posed by the anomalous dispersion of 
the plasmon linewidth in K. The key damping mechanism is shown to be decay into particle- 
hole pairs involving empty states of d-symmetry. While the effect of many-particle 
correlations is small, the correlations built into the "final-state" J-bands play an important, and 
novel, role — which is related to the phase-space complexity associated with these flat bands. 
Our case study of plasmon lifetime in K illustrates the importance of ab initio paradigms for 
the study of excitations in correlated-electron systems. 



Much less is known about how to calculate accurate lifetimes of electronic excitations 
for realistic models of solids than about ground- state and structural properties of these systems. 
Currently a major thrust is focused on the excitations (involving spin, charge, phononic degrees 
of freedom, etc.) in systems of strongly-correlated electrons. The complex band structure 
characteristic of these systems notwithstanding, the prevalent paradigm adopted for their study 
invokes sophisticated treatments of correlation together with minimal representations of the 
band structure [1]. Interestingly, for systems with simpler band structures the puzzles posed by 
a variety of spectroscopic data [2,3] have led to the advocacy of a similar paradigm. This being 
the case, the latter systems offer the possibility — which we pursue in this work — that, on 
account of their relative simplicity, they may be more amenable to ab initio theory. 

As a significant example we have that landmark electron energy-loss spectroscopy 
(EELS) measurements by vom Felde, Sprosser-Prou, and Fink [2,4] revealed remarkable 
anomalies in the excitation spectrum of the alkali metals. For example [5], the dispersion of 
the plasmon linewidth in K (and to a lesser degree in Na) was found to differ qualitatively from 
theoretical predictions obtained on the premise that the effects of the one-electron band 
structure are unimportant, except for the presence of a small gap just above the Fermi surface 
[6]. On this basis it was concluded that (i) the linewidth dispersion is anomalous, and (ii) the 
physics behind the anomaly must be due to strong dynamical short-range correlations [2]. 
Since these correlations vanish in the high-density limit r s — >0, which is far from being 
realized in K, in which the density is low (r s ~ 5), conclusion (ii) seems "natural," and is 
consistent with the absence of single-particle fine structure in the EELS data. 
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In this Letter we show that the anomalous dispersion of the plasmon linewidth in K is, 
in fact, controlled by decay into single particle-hole pairs involving empty states of d- 
symmetry. By comparison with the EELS data we conclude that dynamical many-body 
correlations (which we turn off), do not affect the linewidth dispersion. By contrast, the 
exchange-correlation effects built into the crucial "final-state" J-bands have a profound impact 
on the calculated linewidths. The identification of the physics behind the data is made possible 
by our approaching the problem on the basis of a new all-electron response scheme which is 
framed within time-dependent density-functional theory (TDDFT) [7,8]. The response of 
Kohn-Sham electrons is evaluated for imaginary frequencies, followed by analytic continuation 
to the real axis, which we perform numerically; this procedure leads us to an accurate 
extraction of the natural linewidth. The novel features of our results are related to the phase 
space complexity introduced by the flat J-bands which dominate the physics of the damping; 
this complexity is absent in popular models of correlated-electron systems [1]. 

On first glance, it may seem that the energy scale of the physics involved in plasmon 
damping (~ 4 eV, in K) is too large to be of significance for the current emphasis in the study 
of low-energy excitations and correlations near the Fermi surface. However, in fully self- 
consistent many-body calculations [9] the collective mode feeds back onto the one-particle 
spectrum, and this effect influences, e.g., the value of the quasiparticle renormalization factor at 
the Fermi surface. Moreover, the electronic correlations at low energies are of course affected 
by screening, and this process — usually bypassed via phenomenological on-site interactions — 
is closely related to the physics of the collective mode. 
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In TDDFT the exact linear density-response function ;f (*,*';&>) obeys the integral 
equation [7,8] 

x = x is) +x is) (y+f xc )x , (i) 

where X {s) (x>x / '-> a) ) is the response function for "non-interacting" Kohn-Sham electrons, 
v(x-x') is the Coulomb interaction, and f xc (x,x';Q)) accounts for dynamical exchange- 
correlation effects. For a periodic crystal the Fourier transform of x (s) is given by 



x lk+q,n' e i{ « +6 ' } * k,n) , (2) 



where G is a vector of the reciprocal lattice, n is a band index, the wave vectors k and q are in 
the first Brillouin zone (BZ), and V is the normalization volume. In this representation Eq. (1) 
is turned into a matrix equation which we solve numerically; the plasmon loss corresponds to a 
well-defined peak in Im%g =0 g, =0 (q;co) for a given q . The off-diagonal matrix elements in Eq. 
(2) incorporate the "crystal-local fields" in the response function. This effect turns out to play a 
minor role in the present calculations; thus, it will not be mentioned any further. 

In this work we compute x (s) making use of the local-density approximation (LDA) in 
the evaluation of the exchange-correlation potential V xc (x) [10]. The remarkable success of 
this approximation in the present response context is addressed below. We adopt an all- 
electron approach, which starts from the knowledge of the LDA band structure and wave 
functions [11] in the full-potential linearized augmented plane wave (LAPW) method [12]. 



We aim at disentangling the effects of the one-particle band structure from those of 
dynamical many-body correlations. The TDDFT linear-response framework [8] is well suited 
for this purpose, as it allows us to turn off the latter by setting f xc = — which we do [13]. 
Hence, we concentrate on elucidating the impact of single-particle decay channels. 

Now, an accurate evaluation of the linewidth of a long-lived excitation such as the K 
plasmon poses a non-trivial problem to theory. Indeed, when sampling the BZ in Eq. (2), 
conventional methods introduce numerical broadening via a finite value of J] which typically is 
in the 0.5 eV range — larger than the natural linewidth. We proceed differently. We first 
compute X (s \ an d solve Eq. (1), for imaginary frequencies — this allows us to sample the BZ 
without numerical broadening; J] = in Eq. (2). The physical response % is obtained on the 
real axis — a distance 8 above it — by analytic continuation, which we perform via Pade 
approximants [14]. The power of our method is illustrated in Fig. 1, in which we compare the 
loss function Im^ =0 g, =0 (q;a>) obtained as just outlined [15] with its counterpart obtained via 
the (more standard) direct evaluation of Eq. (2) for real ft)'s (a finite 77 is now required) [16]. 
Figure 1 corresponds to a 16x16x16 wave vector mesh [17]; for the same &-mesh, the real- 
axis approach cannot resolve the linewidth. It is crucial that our calculated loss peak is quite 
insensitive to a decrease of S by an order of magnitude. Since in the upper panel 2d is ~ 1/100 
of the plasmon linewidth, an accurate value of the natural width is thus extracted [18]. 

In Fig. 2 we show a well-converged calculation of the plasmon linewidth dispersion of 
K (solid circles), obtained using a 20x20x20 k -mesh and an energy cutoff of 20 eV in Eq. (2), 
corresponding to the inclusion of ~ 20 valence bands plus core states (the latter states are 
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discussed towards the end of our presentation). Our result is compared with the EELS data of 
vom Felde et al. [2] (empty diamonds). Note that the full- width at half-maximum of the 
plasmon peak, AE l/2 (q) , is given relative to its extrapolated value for q = , a convention also 
adopted by previous authors [2,6]. Clearly, our results are in excellent agreement with 
experiment; since this agreement is obtained for f xc =0, we conclude that the plasmon 
linewidth dispersion of K is not controlled by a dynamical many-body mechanism [19,20]. 

The result of Fig. 2 is striking, as "intuitive" expectations based on the fact that the gap 
just above the Fermi surface at the iV-point is small [21] yield the result shown by the solid line 
in Fig. 2, which corresponds to an evaluation of the dielectric function to second order in an 
empirical pseudopotential [6]. It is apparent that the use of nearly- free-electron states and 
eigenvalues in the evaluation of % {s) breaks down completely. Note that the fact that this 
simple-model result is so far off the EELS data is precisely what led to the proposal of the 
importance of dynamical correlations [2]. 

In Fig. 3 we analyze the role of key "final state" bands entering the evaluation of the 
single-particle response the relevant part of the band structure is shown in the inset of 

Fig. 3 (left panel), in which the shaded strip is the ft) -interval representing all the single-particle 
states which may couple to the plasmon (as determined by the conservation laws of energy and 
crystal momentum) [22]. Now three valence bands (thin solid lines) are needed in order to 
obtain a good dispersion curve for the plasmon energy [23]. Keeping just these bands, our 
plasmon linewidth dispersion curve (solid triangles) agrees well with the result of Ref. [6], 
which is understandable, as the states kept are, for the most part, nearly-free-electron-like. 
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Next, we have the result (solid squares) which incorporates the contribution from three 
additional bands (thick solid lines in the inset). The inclusion of these bands brings about a 
qualitative change in the plasmon linewidth dispersion curve, which is now quite close to the 
EELS data [24] . Clearly, these three bands provide the key decay channels for the plasmon of 
K. It is significant that the bands in question are overwhelmingly of J-character, as evidenced 
by the angular momentum-resolved density of states (DOS) shown on the right panel of Fig. 3. 

The damping channels we have just identified yield a positive plasmon linewidth 
dispersion. We can visualize this point with the aid of the qualitative result that [6] 



where e(q;0)) = 1- v(q)x (s) (q;(o) is the dielectric function (its G = 0,G' = element). The 
first factor in Eq. (3) represents the availability of channels for decay into particle-hole pairs 



disperses upwards [23], a larger portion of the flat bands in the shaded strip becomes available 
to the plasmon as \q\ increases, which naturally leads to an enhancement of the first factor in 
Eq. (3) — thus the positive linewidth dispersion, in agreement with experiment. 

We concluded above that dynamical many-body correlations are unimportant in the 
present problem. However, there is a subtle many-body effect built into our result of Fig. 2; 
the same has to do with the exchange-correlation effects included in the Kohn-Sham response 
X (s) . We emphasize that, although the single-particle states entering Eq. (2) — in particular, 




(3) 



with energy hco p {q) (given by the position of the peak in Imx(q;0))); the second represents a 



polarization effect" from interband transitions taking place at other ft)'s. Then, since CO p (q) 
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the crucial d-bands of interest here — do not have the meaning of quasiparticle states, Eq. (2) is 
rigorous, in the TDDFT linear-response framework [8]; the only approximation we have made 
in the evaluation of x (s) is the LDA. We have performed additional calculations within the 
random-phase approximation (RPA), in which all exchange-correlation effects are left out, 
including those in the band structure — which now corresponds to the Hartree approximation. 
Our RPA results are given on the left panel of Fig. 4, together with the full LDA results. The 
key physical change is that the flat Hartree J-bands are shifted upwards in relation to the LDA 
J-bands of Fig. 3; the former lie almost entirely above the shaded strip in Fig. 4. As a result, 
the plasmon linewidth dispersion curve acquires a much smaller slope, which differs 
significantly from experiment. This finding places our result of Fig. 2 in an even more 
interesting perspective, since it means that the correlations contained in the one-particle-like 
Kohn-Sham states of d-symmetry play a non-trivial role — via % {s) — in the explanation of the 
plasmon damping mechanism. These correlations are usually not dealt with explicitly in many- 
body models of correlated electrons. 

Finally, we note that K has a polarizable core, which impacts the damping process via 
the second factor in Eq. (3). The "local-orbital" extension of the LAPW method [12] allows us 
to treat the contribution to x (s) from higher-lying core states, mainly derived from the 3p 6 
atomic states (whose threshold is at 16 eV), on the same footing with that from the valence 
states [25]. The circles (squares) on the right panel of Fig. 4 are the result obtained in the 
presence (absence) of this contribution to It is apparent that the absolute value of the 

linewidth is affected markedly by the availability of the core excitations. In fact, their effect on 
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the slope of the linewidth dispersion curve is even larger than that of the many-particle 
correlations [23]. Of course, this effect was included in our central result of Fig. 2. 

In summary, we have explained the "anomalous" dispersion of the linewidth of the K 
plasmon via ab initio TDDFT-based calculations. The key mechanism was shown to be decay 
into particle-hole pairs involving empty states of d-symmetry. Our results highlight the role of 
the exchange-correlation effects built into the LDA single-particle band structure — these 
effects become important because the key final-state bands are flat. It is the phase-space 
complexity associated with these bands which renders simple models inapplicable. Ab initio 
paradigms should also play an important role in the study of excitations in materials with 
complicated band structures near the Fermi surface [26] . 
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FIGURE CAPTIONS 



Fig. 1. Loss function for K for q = (1,1,0)(2tt / (16<3 )), where a =5.23A is the lattice 
constant. Thick solid line: Im^ obtained from Eqs. (1) and (2) for imaginary ft)'s, 
followed by analytic continuation to a distance 5 above the real- ft) axis [15]. Thin 
solid line: Im^ obtained directly for real ft) 's, using a broadening parameter r\ in 
Eq. (2). Top (bottom) panel corresponds to 8 = r\ = lmeV (lOmeV ). 

Fig. 2. Plasmon linewidth dispersion for K. Comparison of our theoretical results with the 
EELS data of Ref. [2] (diamonds), and the theoretical results of Ref. [6] (solid line). 
Theory is for (1,1,0) propagation; the EELS data are for polycrystalline K. 

Fig. 3. Left panel: Plasmon linewidth dispersion obtained upon keeping 3 (triangles) and 6 
(squares) valence bands in x [24]. Inset: LDA band structure of K; the arrow 
indicates the value of ft) p (0)(see text and [22]). Right panel: Calculated DOS for K 
— total DOS and contributions from states of s,p, and d- symmetry [11]; the zero of 
energy is the Fermi level. 

Fig. 4. Left panel: Plasmon linewidth dispersion obtained on the basis of RPA (diamonds) 
and LDA (circles) single-particle response x (s) and respective band structures. The 
inset shows Hartree bands used in the RPA case; note the impact of exchange- 
correlation effects on the location of the J-bands (cf. Fig. 3). Right panel: 
Calculated plasmon linewidth dispersion in the presence (circles) or absence 
(squares) of the core contribution to the Kohn-Sham response x (s) ■ See text. 
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